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In practical implementations of density-functional theory, the only term where an orbital description 
is needed is the kinetic one. Even this term in principle depends on the density only, but its 
explicit form is unknown. We provide a novel solution of the iV-representability problem for an 
extended system, which implies an explicit form for the Kohn-Sham kinetic energy in terms of the 
density. Our approach is based on a periodic coordinate mapping, uniquely defined by the Fourier 
coefficients of the metric. The density functional is thus expressed as an explicit functional of the 
metric tensor: since ./V-representability is enforced, our constructive recipe provides a variational 
approximation. Furthermore, we show that our geometric viewpoint is quite naturally related to 
the electron localization function (ELF), which provides a very informative analysis of the electron 
distribution. Studies of ELF, as obtained from accurate Kohn-Sham orbitals in real materials, 
allow an appraisal of the variational approximate density functional. We show that the value of 
an approximate functional — either the present geometric-based one or some previous ones based 
on different constructive recipes — strongly depends on the nature of the chemical bonding in the 
material. 
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I. INTRODUCTION 

The |Gelebrated basic tenet of density-functional theory 
(DFT)ETa states that an exact description of a many- 
electron system is in principle possible in terms of a 
single scalar field, namely the electron density n(r). 
The iV-electron wavefunction contains instead redun- 
dant information; in extendend systems, it does not even 
have a well defined thermodynamic limit. However, the 
Hohenberg-Kohn theorernEl (upon which DFT is based) 
does not provide a constructive scheme: for any given N, 
the exact functional is indeed accessible only through the 
many-body wavefunction. 

The enormous success of DFT resides in approx- 
imate schemes which are constructive, and do not 
make explicit recourse to the many-body wavefunc- 
tion. These schemes — implemented using first-principles 
ingredients — have proved over the years their astonish- 
ingly accurate predictive power for many physical proper- 
ties, in many different materials. In all these schemes the 
only wavefunction needed is a wavefunction of noninter- 
acting electrons, which is uniquely defined by the mani- 
fold o£-jtl^3 occupied Kohn-Sham (KS) single-particle or- 
bitals,Ba or equivalently by the KS one-body reduced 
density matrix. The eponymous density functional F[n], 
Eq. (nOl) below, is the sum of a few terms: all of them 
but one are almost invariably approximated as explicit 
functionals of the density. The only term where the KS 
orbitals (or the density matrix) are actually needed is T s , 
the kinetic energy of the noninteracting system, which is 
a functional of the density in an implicit way. This qual- 
itative difference is of course responsible for most of the 
computer workload in practical calculations, hindering 
amongst other things the linear scaling of computations 
with the size of the system. There is therefore a quest for 



approximate (though accurate enough) expressions for T s 
as explicit functionals of the density. 

Historically, the first approximate form for the ki- 
netic energy of a system of noninteracting electrons i: 
terms of their density is the Thomas-Fermi (TF) onej 
which predates DFT by several decades and is a very 
crude one; it is however exact for an extended system 
of free electrons. We focus here on a different class 
of approximations, which are — at variance with TF — 
variational: this feature is intimately linked to the prob- 
lem of iV-representability. In fact, whenever the approx- 
imate T s coincides with the true kinetic energy of an 
arbitrary independent-electron wavefunction, then the 
variational theorem ensures that it must be no smaller 
that the exact T s of the given system. The requisite of 
iV-representability is equivalent to requiring that T s ob- 
tains from a density matrix which is idempotent. There 
is clearly an infinity of idempotent density matrices, all 
yielding the same given density: amongst these infinite 
solutions of the ./V-representability problem, one searches 
for the one having the lowest T s at the given density. 

In order to provide an explicit approximate (and vari- 
ational) expression for T s one has, first of all, to pro- 
vide a constructive recipe which, starting from a given 
density, produces an idempotent density matrix. Ex- 
plicit solutions of this problem have been provided by 



several authors i 
quote Harrimaj 
and coworkers 



I the literature, amongst whom we only 
^umbach and Maschke,0 and Ludena 



The present paper may thus be con- 
sidered as a continuation and a generalization of this ear- 
lier work, where the elements of novelty are basically the 
following, (i) At variance with previous work, we are in- 
terested in extended systems. We therefore solve the N- 
representability problem for a system of N noninteracting 
electrons in a box of volume V, where periodic (Born- von 



Karman) conditions are assumed at the boundary. Our 
approximate solution coincides with the exact one for the 
electron gas. (ii) Our construction uses in an essential 
way a periodic coordinate m appi ng, in the same spirit 
as the one advocated by GygiliilO in electronic structure 
calculations. This provides an elegant and symmetric 
treatment: some of the results obtained in the previous 
literature assume here an interesting geometric meaning. 
The electronic energy is variationally expressed in terms 
of the metric tensor as the independent variable, (iii) 
Our geometric approach naturally partitions T s into the 
sum of two terms: very roughly speaking "bosonic" and 
"Pauli" . We show that these two terms coincide with the 
volume integrals of two local functions which are used in 
the literature as the. main ingredients of the electron lo- 
calization functionl3'li-3 (ELF). 

Our solution of the ./V-representability problem — as 
well as the explicit approximate density functional based 
upon such solution — is therefore fundamentally linked in 
a very natural way to the ELF concept. Since ELF pro- 
vides a very informative analysis of the electron-distri- 
bution, the published results for real materialsEjta — 
obtained a-posteriori from accurate density matrices — 
help understanding what is good and what is bad in the 
approximate forms of T s . We find that the quality of 
the approximation provided by our constructive recipe 
strongly depends on the kind of bonding involved in the 
many-electron system. We give evidence that both the 
Zumbach-Maschke recipdJ and our own one provide a 
reasonable T s for simple metals in the pseudopotential 
approximation, while they significantly overestimate T s 
whenever covalent bonding is present. The possible di- 
rections for improvement on this point are sketched. 

The present paper is organized as follows. In Sec. II we 
present our novel solution to the iV-representability prob- 
lem, essentially based on a coordinate mapping, which 
transforms a reference uniform system into the actual 
nonuniform one. In Sec. Ill we apply such solution to 
the construction of a functional which is a variational 
approximation to the exact one: this functional is an 
explicit functional of the metric tensor. In Sec. IV we 
review the foundamental properties of ELF, a powerful 
tool used in the quantum-chemistry community to ana- 
lyze the electron distribution in various real systems: we 
show that ELF is intimately related to basic features of 
the functional. Finally, in Sec. V we outline our conclu- 
sions and perspectives, based on some ELF analyses for 
real materials. 



II. CURVILINEAR COORDINATES AND 
7V-REPRESENTABILITY 

For a system of independent electrons in a closed-shell 
configuration the wavefunction is a single determinant: 
knowledge of the one-particle reduced density matrix is 
equivalent to a complete knowledge of the wavefunction. 



The spin-integrated matrix p(r,r') is twice a projector, 
which indeed projects over the (doubly occupied) one- 
particle orbitals. We consider a system of N electrons in a 
box of volume V, obeying periodic boundary conditions. 
The average density is no = N/V, the density is n(r) = 
p(r, r), and the idempotency condition is written: 



[dv"p(r,r")p(v",v') = 2p(x 

Jv 



rO- 



ll) 



We start with a homogeneous system of N non- 
interacting electrons at the same density no, for which we 
use £ as a space coordinate. For this system the canoni- 
cal orbitals are, by symmetry, the plane waves e lk '^/VV, 
where k; are the reciprocal vectors determined by the 
boundary conditions. By choosing to occupy the N/2 
orbitals of lowest energy, the density matrix is: 



N/2 



po(^o = 4E eik,( ^' ) 



(2) 



1=1 



which is obviously idempotent and yields a constant den- 
sity. In the thermodynamic limit (N — ► oo and V — > oo 
at constant no) the k; set becomes dense. Occupying 
the k-vectors within the Fermi sphere (|k| < fep) Eq. (g) 
yields the well known electron-gas result: 



po(€,C) 



n - 



3 ji(fc F lg-g / l) 



(3) 



At this point we introduce a generic curvilinear co- 
ordinate precisely of the same kind as introduced-Jpji 
Gygi in the field of electronic-structure calculations lM3 
We therefore define a twice differentiable invertible map 
£ — > r(£), periodic over V, whose Riemannian metric 
tensor is: 



m 



dr k dr k 
WW' 



(4) 



Summation over repeated indices is understood through- 
out. A generic plane wave of momentum k is transformed 

as: 



— e ik * 



Xk(r) = -==g~ 



■(r)e 



;k£( r 



(5) 



where g = det{gij}, and g~? is the Jacobian \d£/dr\ 
of the inverse transformation. Notice that the orbitals 
Xk( r ) have a k-independent density, and are therefore 
"equidensity orbitals" in Harriman'sEl nomenclature. The 
density matrix in the new coordinates is: 



p(r,r')=A)(£(rU(r')) = 



N/2 

-g-i(r)g-i{r')'2^e- 
i=i 



r /)y^ e ik,[C(r)-$(r')] > 

1=1 

The corresponding transformed density is 



(6) 



n(r) = n Q g »(r). 



(7) 



In the novel coordinates we thus have a nonhomogeneous 
system, with the same average density as the homoge- 
neous one, and whose density matrix is idempotent by 
construction. 

We are now ready to attribute physical content to the 
above mathemathics. Suppose that the density n(r) of 
an electronic system is given. Then we may look for a 
coordinate transformation £ — ► r(£) which maps the 
uniform density into the given density: a necessary and 
sufficient condition is Eq. (ffl). The solution is nonuniquc, 
since several different maps share the same Jacobian g~?: 
we will argue below about an optimal solution, usipffj-a 
variationally adaptive metric in the sense of Gygi.O : E£l 
Replacement into Eq. (0) yields the explicit form: 



p(r,r') 



N 



n 2 (r)n 2 (r') 



N/2 

E< 



iki[£(r)-£(y)] 



(8) 



In one dimension the solution of Eq. ([7]) is unique, 
and we get here the periodic analogue of the Harriman 
construction.!] In three dimensions, our-jesult is related 
to the work of Zumbach and MaschkeJ3 the differences 
being that we deal with periodic systems, and we pro- 
vide a more general explicit construction. A coordinate 
mapping, similar in spirit to the present one (and called 
"local scaling transformation"), has been previously in-i 
troduced by Ludeha and coworkers for spherical atoms.cl 



III. DENSITY FUNCTIONAL 

The energy of the electronic system k^the external 
potential v cx t is written, within DFT, asi 



E[n] = / drn(r)v ext (r) + F[n}; 
Jv 

m r i 1 f , f , / nMn(r') 
F[n] = T s [n] + - dr dr' , V ; \.> 



(9) 

£«M (io) 



where atomic Hartree units have been used. As already 
anticipated, basically all the available constructive ap- 
proximations to DFT provide E xc as an explicit func- 
tional of the density n(r), while instead the kinetic energy 
term XL is: 



T s = ~ /drV r Vw>(r,r'; 



(11) 



The ground electronic energy, Eq. (0), is therefore an 
explicit functional of the density matrix, which has to be 
minimized under the constraints of idempotency, Eq. ffl), 
and electron number. 

Replacement of our ansatz density matrix, Eq. (ph, in 
the above expressions provides an upper bound to the 
electronic energy, explicitly expressed solely in terms of 



the density and of the metric. The approximate kinetic 
energy is: 



JV/2 

t = J2 *|V Xki (r) 



(12) 



Using then Eq. (7) of Ref. |l2j, the expectation value of 
the kinetic energy over a Xk orbital is the sum of two 
positive terms: 

(Xk|T| Xk ) = ^ J v d£g ij + ± J^i A^Aj, (13) 



where the "gauge potential" is: 

1 d In q 

A . — _ 

1 ~ 4 d? ' 



(14) 



After summing over the N/2 doubly occupied states, we 
get 



f 8 [n] = f P [n]+T B [n], 



(15) 



where the reason for the notations will be clear in a mo- 
ment. 

Using Eq. (0), we cast the gauge term as: 



r B [n] = ^jdrg-i(r) |V 5 (r)| 2 = \ j dv 



|Vn(r) 
n(r) 



(16) 



In the latter expression, we notice that the metric for- 
mally disappeares from the gauge term, which is in- 
deed identical to the so-called von Weizsacker energy 
functional!] This energy coincides with the kinetic en- 
ergy of a system of noninteracting bosons in their ground 
state, having the given density n(r): with this specific 
meaning, we may refer to Tb as to the "bosonic" energy. 
It is easy to prove that Tb is a lower bound to the kinetic 
energy T s of a system of noninteracting fermionsp and 
coincides with T s only in the trivial case N = 2, where 
the Pauli principle has no effect (we are considering sin- 
glets only). In all the interesting cases, there is an excess 
kinetic energy Tp due to the Pauli principle: 



T P = l [ drV r V r ,p(r,r') - 1 / dr 

Z JV r=r' ° JV 



|Vn(r)| 5 
n(r) 



(17) 



Our expression for Tp [n] as a function of the metric shall 
therefore be a variational approximation to the true ex- 
cess Pauli energy Tp. 

The sum over the occupied states in Tp is most easily 
evaluated if we assume a cubic box. If we define Eq as 
the kinetic energy of the homogeneous system: 



E 



N/2 

Ei k <i 2 > 



(18) 



it is then easy to recast Tp [n] as: 



f P [n] = 



-£•0 
3V~ 



drg 2tr{ 5 y }. 



(19) 



We further notice that in the thermodynamic limit one 
has: 



10 



E = —N4 = c F Vnl. 

(3^ 2 )f. 



(20) 



where cp 

Putting all the previous formulas together and approx- 
imating T s in Eq. (nfl) with T s , we obtain an approximate 
F[n] as an explicit functional of the density and of the 
metric tensor. Since the density — owing to Eq. ffik — is 
in turn a function of the metric tensor, we use the latter 
as the independent variable. Eventually, the electronic 
energy of the system, Eq. (0), is a variational explicit 
functional of the metric tensor g u (r). This functional 
can be regarded as the_periodic analogue of the one of 
Zumbach and Maschke,0 expressed in more compact form 
in terms of a different variable. Furthermore the explicit 
occurrence of the periodic metric in the Tp term makes 
feasible an adaptive optimization n£ the metric. Upon 
closely following Gygi's approach,Ej the Fourier coeffi- 
cients of the periodic metric are the natural variational 
parameters of the problem. 

Finally, we end this section just noting that in the 
trivial case u cx t(r) = all of the kinetic energy is due 
to Tp, since the density is constant and Xb vanishes. 
Furthermore the metric is the identity and the approxi- 
mate kinetic energy equals the exact one: T s = T s = Eq, 
Eq. (p0|). The approximate functional F[n] coincides with 
the exact one, including itS|-exchange-correlation term if 
the exact electron-gas dataO are used therein (as usual). 
This suggests that the approximate functional should 
work reasonably well for a system close enough to the 
electron gas, such as a simple metal within a pseudopo- 
tential scheme. 



IV. ELECTRON LOCALIZATION FUNCTION 

The kinetic energy T s can be thought of as the integral 
over V of a kinetic energy density r(r). It is well known 
that the expression for t(t) is not unique: we use the 
form suggested by Eq. ([fl]), namely, 



r(r)= -V r V r -p(r,r') 



(21) 



which is everywhere positive, and -coincides with the 
choice made in the ELF literatures By analogy, one 
defines the Pauli excess energy density, after Eq. (17), 
as: 



rp(r)= iv r V r ,p(r,r'; 



1 |Vn(r)| 2 



(22) 



This function is the main ingredient of ELF ,113 in the 
formulation due to Savin et al.S-3 who write the function 



£(r) 



cp ni(r) 



(23) 



This function by design takes values between zero and 
one: several of its features are remarkable. In the homo- 
geneous electron gas, owing to Eq. (p0|), the ELF equals 
1/2 at any density. In a nonhomogenous system £ (r) as- 
sumes values close to its upper bound 1 in the regions 
of space where there is an high probability of finding a 
pair of elections with antiparallel spins (or an isolated 
electron) :liijlia with this meaning, we may say that £ (r) 
close to 1 characterizes space regions where the electron 
distribution is "bosonipli. Conversely, £(r) is close to 
in low-density regions EB 





0.0 0.5 1.0 

FIG. 1. Pseudocharge density contour plot (upper panel) 
and the corresponding ELF (lower panel) for bulk silicon in 
the [110] cristalline plane. The grey-scale is also shown: dark 
regions correspond to large ELF values. The maximum ELF 
value at the bond center is 0.96. 





0.0 0.5 1.0 

FIG. 2. Pseudocharge density contour plot (upper panel) 
and the corresponding ELF (lower panel) for bulk aluminum 
in the [100] cristalline plane. The maximum ELF value in the 
interatomic region is 0.61. 

Detailed ELF analyses have-heen performed for many 
systems of chemical interest .li3~EZI For heavy atoms, the 
ELF perspicuously localizes in space the different elec- 
tronic shells: furthermore, when the valence shell is con- 
sidered, the ELF provides a very meaningful description 
of the chemical bond for many classes of compounds.li.lt3 

For the sake of simplicity, we consider only the class 
of sp-bonded materials which are accurately described 



within a norm-conserving pseudopotential schemes this 
class includes simple metals, covalent semiconductors, 
simple ionic solids, and many other disparate materials. 
We have in this case by construction only a single shell 
(the valence sp one of each atom involved) , and the ELF 
allows a very meaningful analysis of the chemical bond- 
ing. Perhaps the most spectacular performance is the 
ability to perspicuously distinguish in a very clearcut way 
between metallic bonding and covalent bonding. This is 
illustrated in the lower panels of Figs. (|]) and (g), where 
we plot the function £ (r) for two paradigmatic crystalline 
materials: respectively, silicon and aluminum. 

In the covalently bonded system of Fig. |l], the bond 
electron pairs forming the typical "zig-zag" chain in the 
[110] plane are clearly visible. The very dark regions in 
the lower panel indicate in fact the strong bosonic char- 
acter of the charge density in the bond region. Actually, 
£(r) attains the maximum value of 0.96 at the bond cen- 
ter, thus indicating an extremely strong pairing between 
opposite-spin electrons. 

A completely different picture emerges for our sim- 
ple metallic system, Fig 0. The ELF plot in the lower 
panel shows — outside the core regions — a large grey area, 
which correspond to a jellium-like (or Thomas-Fermi) 
ELF value. Actually, the maximum value attained by 
£ (r) is only 0.61. This is agreement with the usual pic- 
ture of a simple metal, where the (pseudo) electrons are 
excluded from the core region, while behaving essentially 
as free particles in the rest of the material.E3 

Comparison of the two ELF plots provides therefore 
the most significant and perspicuous visualization of 
the important qualitative difference between the covalent 
bond and the metallic one. In the two classes of mate- 
rials the Pauli principle plays quite a different role. At 
variance with the ELF, the corresponding charge density 
plots, (upper panels in Figs [fl and 0) are much less in- 
formative, and do not qualitatively discriminate between 
the two different kinds of chemical bonds. 



V. CONCLUSIONS AND PERSPECTIVES 

The experience gained in investigating the ELF in real 
materials helps understanding the meaning and the lim- 
its of the approximate density functional such as the 
Zumbach-Maschkeu one, as well as of the generalization 
proposed here in Sec. III. 

Our explicit ansatz of Eq. (11) leads to the Pauli excess 
energy Tp of Eq. (|l^). It is interesting to see the con- 
sequences for the ELF, since the ansatz clearly leads to 
replacing the Pauli excess energy density, Eq. (E2h, with 



f P (r) = ^ g-h tr{gV] = ippnj <T* tr{g»}, (24) 



where the thermodynamic limit, Eq. (|20|), has been used. 
Considering now the inequality 



\trW j } > [det{ 5 «}]i = g-i = ^R (25) 

we get for the approximate Pauli energy density the lower 
bound: 



■?p( r ) > cp n 3 (r). 



(26) 



Comparing with the ELF definition, Eq. (p3h, one easily 
realizes that even the optimal choice of the metric tensor 
will unavoidably provide values of £ (r) which are smaller 
than 1/2 everywhere. 

In this work we have shown that a fundamental link 
exists between ./V-representability, approximate explicit 
functionals, and ELF. This link emerges very naturally 
within the geometric approach upon which our work is 
based. The same geometric approach, however, also in- 
dicates very clearly the limits of the approximate form 
of the kinetic energy for real materials which we have 
found here. In fact the final considerations of the pre- 
vious paragraph imply that our constructive recipe, as 
well as the previous one of Zumbach and Maschke,Q are a 
good approximation only for systems where the bonding 
is metallic, while it necessarily overestimates the kinetic 
energy (and the total energy) whenever covalent bond- 
ing is present. Looking more closely, this major limita- 
tion owes to the occurrence of equidensity orbitals in our 
ansatz density matrix, Eq. (|6j), which occurrence can be 
further traced back to the choice of the uniform electron 
gas as the reference system upon which we perform the 
coordinate transformation. 

This naturally suggests the directions for improve- 
ments: one should start from a reference model sys- 
tem other than the uniform electron gas, having instead 
some covalent bonding features already built in. Inter- 
estingly, the use of a model reference system designed 
to reproduce — after coordinate mapping — some desirable 
features of the real one has been pronosed in the most 
recent work of Ludeha and coworkers. B'Ej These authors, 
however, focus on a spherical atom having several elec- 
tronic shells: here instead we are discussing a condensed 
system with only one valence shell, within a pseudopo- 
tential scheme. 
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